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Abstract 

This  paper  presents  a  new  look  at  Weispfenning’s  method  of  quantifier 
elimination  by  virtual  term  substitution  and  provides  two  important  im¬ 
provements.  Virtual  term  substitution  eliminates  a  quantified  variable  by 
substituting  formulas  in  the  remaining  variables  for  each  atomic  formula 
in  which  the  quantified  variable  appears.  This  paper  investigates  the  poly¬ 
nomials  that  arise  in  substitution  formulas  Weispfenning  proposed  and, 
based  on  this  examination,  provides  a  simpler  substitution  for  the  general 
case,  and  alternate  substitutions  for  several  commonly  occurring  situtions. 
Providing  alternate  substitutions  allows  virtual  term  substitution  to  make 
choices  that  produce  simpler  output. 


1  Introduction 


Quantifier  elimination  for  elementary  real  algebra  is  a  fundamental  problem 
in  symbolic  computing.  The  great  potential  utility  of  quantifier  elimination 
algorithms  is,  however,  offset  by  an  equally  great  theoretical  and  practical  com¬ 
plexity.  Thus,  the  search  for  improved  algorithms  that  are  capable  of  solving 
interesting  problems  in  a  reasonable  amount  of  time  and  space  is  important. 

One  successful  result  of  this  ongoing  search  is  Weispfenning’s  method  of  quan¬ 
tifier  elimination  by  virtual  term  substitution  [4,  6,  5].  This  method,  restricted 
to  formulas  that  are  linear  or  quadratic  in  the  quantified  variable,  has  been 
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implemented  in  the  Redlog  system  [2],  and  has  been  applied  successfully  to  a 
number  of  practical  problems  (e.g.  [7]). 

The  strength  of  quantifier  elimination  by  virtual  term  substitution  is  that  its 
complexity  is  relatively  unaffected  by  the  number  of  parameters  -  i.e.  un¬ 
quantified  variables  —  in  the  problem.  The  method  has  two  weaknesses:  First, 
applying  the  method  iteratively  to  eliminate  several  quantified  variables  may 
not  be  possible  because  eliminating  one  variable  may  increase  the  degrees  of  re¬ 
maining  variables,  thus  violating  the  degree  restrictions.  Second,  quantifier-free 
equivalent  formulas  produced  by  the  method  tend  to  be  very  large,  even  when 
simple  quantifier-free  equivalents  exist. 

The  purpose  of  this  paper  is  to  provide  a  new  perspective  on  virtual  term  sub¬ 
stitution,  and  to  apply  this  new  perspective  to  help  address  the  method’s  weak¬ 
nesses. 


1.1  What’s  new 

Virtual  term  substitution  is  based  on  rewriting.  A  formula  3a;[F’]  is  transformed 
into  an  equivalent  formula  in  which  x  does  not  appear  by  combining  many  copies 
of  F  in  which  atomic  formulas  are  substituted  by  more  complex  formulas.  In  this 
paper  we  give  a  new  analysis  of  which  polynomials  appear  in  these  substitution 
formulas  and  why.  Using  this  analysis,  we  provide  an  improved  substitution 
(Section  3)  for  the  algorithm’s  general  case:  improved  in  that  fewer  atomic 
formula  and  fewer  distinct  polynomials  appear  in  the  output  formula.  Generi- 
cally,  the  improved  method  produces  a  formula  in  which  the  set  of  polynomials 
occurring  is  a  proper  subset  of  those  occurring  in  the  original  method’s  output. 

Also  using  our  initial  analysis,  we  provide  alternate  substitutions  which  may 
be  used  in  certain  situations  —  allowing  the  method  to  evaluate  alternatives 
and  choose  the  best  substitution  (see  Section  4).  In  particular,  substitutions 
which  use  lower  degree  polynomials  or  which  allow  obvious  simplifications  can 
be  chosen. 

Two  detailed  examples  applying  both  improvements  are  given  in  Section  5.  The 
ideas  presented  here  have  not  yet  been  implemented,  so  experimental  data  on  a 
wide  range  of  benchmark  problems  cannot  be  given. 

Finally,  Section  6  applies  the  characterization  from  previous  sections  of  the 
polynomials  appearing  in  substitutions  to  determine  a  new  bound  on  the  degree 
of  any  irreducible  factor  of  a  polynomial  appearing  in  the  formula  resulting  from 
eliminating  a  variable. 
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2  Evaluating  a  formula  at  a  quadratic’s  roots 


Suppose  we  want  to  eliminate  x  from  the  formula  3 x[f  =  0  A  F],  where  /  = 
ax2  +  bx  +  c  and  each  atomic  formula  in  F  is  of  the  form  g  a  0.  Theorem  2.1 
of  [6]  describes  the  virtual  term  substitution  approach  to  solving  this  problem. 
Virtual  term  substitution  starts  by  substituting  (— b±  \/D)/2a  for  x  in  F.  Each 
atomic  formula  g{{— b±  \/D)/2a)  a  0  is  then  replaced  by  an  equivalent  formula 
without  radicals. 

For  the  rest  of  this  section  we  let  /  =  ax2  +  bx  +  c,  D  =  b2  —  4 ac,  a+\  = 
{— b  +  '/D)/2a,  a_i  =  (— b  —  \/D)/2a ,  and  let  g  be  an  integral  polynomial 
of  positive  degree  n  in  x.  First  we  will  prove  some  simple  results  concerning 
g(a±i).  Then  we  will  restate  the  substitution  described  in  Theorem  2.1  of  [6]. 


2.1  Evaluating  g  at  the  roots  of  / 

The  following  results  provide  several  characterizations  of  the  evaluation  of  a 
polynomial  at  a  root  of  a  quadratic.  They  are  the  basis  of  this  paper’s  discussion 
of  virtual  term  substitution. 


Lemma  1  For  k  >  1,  (— b  +  y/b2  —  4 ac)k  =  2fc_1  ( Uk  +  Vk'Jb2  —  4ac),  where 
Uk  and  Vk  are  integral  homogeneous  polynomials  of  total  degree  k  and  k  —  1 
respectively,  such  that  Uk  —  bVk  has  even  integer  content. 


PROOF.  We  proceed  by  induction  on  k.  The  lemma  clearly  holds  for  k  =  1 
since  U\  =  — b ,  V\  =  1,  and  U\  —  bV i  =  —2b.  Suppose  the  lemma  holds  for  some 
k.  Then 

(-b  +  Vb2  -  4ac)k+1  =  2k~1(Uk  +  Vk Vb2  -  4 ac)(-b+Vb2  -  4ac) 

=  2fc"1  ((-b(Uk  -  bVk)  -  4acl4)  +  {Uk  -  bVk)\/b 2  -  4ac) 

By  supposition,  Uk  —  bVk  =  2 A  for  some  integral  polynomial  A  of  total  degree 
k ,  which  clearly  must  be  homogeneous.  Therefore, 

— b{Uk  —  bVk)  —  4acI4  =  —  b(2A)  —  4acl4  =  2{—bA  —  4acl4) 

which  has  total  degree  k  +  1  and  is  also  clearly  homogeneous.  So  {—b  + 

\/b 2  —  4ac)fe+1  =  2kUk+i  +  2kVk+i\fSr~-4ac,  where  Uk+i  =  —bA  —  4acVfc  and 
Vfc+i  =  A.  Finally,  we  note  that 

Uk+ 1  —  bVk+i  =  —bA  —  4acf4  —  bA=  —2  bA  —  4acVk 

which  clearly  has  even  integer  content.  [H 
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Theorem  1  Formal  substitution  of  x  =  ap  into  g{ x)  yields  (a*  +pb* \/D) / (2a™) , 
where  a*  and  b*  are  integral  polynomials. 


Theorem  2  Let  a*  and  b*  be  as  above  and  let  5  =  n  mod  2.  res x{f,g)  =  (a*2  — 
b*2D)/ (4a™)  and  sgn(g(a+i)g(a-i))  =  sgn(a*2  -  b*2D)  =  sgn{asresx(f,  g)). 


PROOF.  By  a  well-known  theorem  (see  for  example  [1])  res x{f,  g)  =  anWlLlg{ai), 
where  m  =  degx(f),  n  =  degx(g)  and  the  roots  of  /  are  a\, . . . ,  am.  So: 

res  x(f,g)  =  ang(=^)g(=^) 

_  r.n  (  a*  -\-b*  y/~D  \  (  a*  —b*  y/D  ^ 

~  a  y  2 an  J  \  2 an  ) 

=  (a*2  —  b*2D)/(Aan) 

The  rest  follows  easily.  Q 


Theorem  3  Let  prem(c/,  /)  =  rx  +  s.  Then 

1.  b*  =  r  =  psci(g,  f), 

2.  a*  =  2 as  —  br  =  —rfx{—s/r),  and 

3.  4a”resx(/, g)  =  4 a(as2  —  srb  +  r2c)  =  a*2  —  b*2D. 


PROOF.  Since  prem(g,  /)  =  rx  +  s,  we  have  an  1g  =  Qf  +  rx  +  s,  where  Q  is 
a  polynomial.  So 


an-lg{^D) 

an-lg{=b±^D) 

9(^) 


+  S 


Since  prem(g, /)  =  sresi (g,f),  r  =  psc That  2 as  —  br  =  —res X(rx  + 
s ,  fx)  =  —rfx(—s/r)  and  4a"resx(/,  g)  =  a*2  —  b*2D  =  4a(as2  —  srb  +  r2c)  are 
easily  checked  by  simple  calculations.  F] 


2.2  Virtual  term  substitution  with  a  quadratic  constraint 

We  are  now  ready  to  restate  the  substitutions  from  Theorem  2.1  of  [6],  which 
eliminate  x  from  a  formula  of  the  form  3x[f  =  0  A  F], 
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Theorem  4  (Weispfenning)  Let  f  =  ax2  +  bx  +  c  and  let  g  be  an  integral 
polynomial  of  degree  n  in  x.  Let  6  =  n  mod  2  and  let  R  =  4 a(as2  —  srb  +  r2c). 
Assuming  a  ^  0  A  b2  —  4 ac  >  0, 

1.  g(ap)  =  0  •<=>■  pra*  <  0  A  R  =  0. 

2.  g(ap)  ^  0  pra*  >  0  V  R  ^  0. 

3.  g(ap)  <  0  ■<==>  a*  a5  <  0  A  R  >  0  V  pras  <  0  A  (a* a5  <  0  V  R  <  0) 

4-  g{ap)  <  0  a*  a5  <  0  A  R  >  0  V  pra 5  <  0  A  R  <  0 


Corollary  1  Given  the  assumptions  above, 

g{otp)  <  0  a*as  <  0  A  R  >  0  V  pra5  <  0  A  {a* a5  <  0  V  R  <  0). 

PROOF.  Note  that  3)  from  above  can  be  written  equivalently  as  a* a5  <  0A R> 
0  V  pra5  <  0  A  ( a*a 5  <  0  A  R  =  0  V  R  <  0).  Recall  that  R  =  a*2  —  r2D,  so 
R  <  0  =>  r/  0  and  R  =  0  A  r  =  0  =>  a*  =  0.  So,  pra 5  =  0  is  inconsistent 
with  both  a* a5  <  0  A  R  =  0  and  R  <  0.  Q 


There  is  some  reason  to  consider  this  alternative  substitution  for  g  <  0,  since 
a  simplifier  that  removes  inconsistent  subformulas  might  more  easily  recognise 
an  unsatisfiable  branch  with  pra 5  <  0  rather  than  pra5  <  0  -  if  r  was  a  sum 
of  squares,  for  instance. 

Let  Fap  be  the  formula  obtained  by  replacing  each  atomic  formula  g  a  0  with 
the  appropriate  formula  from  Theorem  4.  Let  F_c/b  be  the  formula  obtained 
by  replacing  atomic  formula  g  a  0  with  res(bx  +  c,g),  noting  that  if  b  ^  0, 
res(bx  +  c,g)  =  b2g(—c/b). 

Theorem  5  (Weispfenning)  Under  the  assumption  a^OV&yfOVcy^O 
3 x[f  =  0  A  F]  a  =  0  A  b  ^  0  A  F_c/b  V  a  yf  0  A  b2  —  Aac  >  0  A  (Fa_1  V  Fa+1). 

This  gives  us  quantifier  elimination  for  formulas  of  the  form  3 x[f  =  OAF]  under 
the  assumption  that  a,  b  and  c  do  not  vanish  simultaneously. 


3  Evaluating  a  formula  near  a  quadratic’s  roots 


Theorem  3.1  of  [6]  gives  a  method  for  eliminating  x  from  a  formula  3a;F  that 
does  not  necessarily  have  an  equational  constraint,  provided  that  all  irreducible 
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factors  of  polynomials  in  F  (recall  that  atomic  formulas  are  normalized  to  g  a  0) 
have  degree  at  most  two.  The  approach  is  based  on  introducing  the  positive 
infinitesimal  e  in  the  substituted  expressions  and  formal  substitution  of  —  oo. 

The  idea  is  as  follows,  there  is  an  x  satisfying  F  if  and  only  if  F  is  satisfied  at 
i  =  a  or  i  =  a  +  f  for  some  real  root  a  of  a  polynomial  in  F,  or  at  x  =  —  oo. 
This  is  clear  because  the  truth  value  of  F  can  only  change  as  x  passes  through  a 
root  of  the  left-hand  side  of  some  atomic  formula.  Thus,  3x[F]  is  equivalent  to 
the  disjunction  of  F  evaluated  at  each  of  these  candidate  points.  Weispfenning 
improves  on  this  by  showing  that  if  a  polynomial  /  only  appears  in  the  atomic 
formulas  /  =  0  or  /  <  0,  the  point  x  =  a  +  e  where  /(a)  =  0  does  not  need  to 
be  tested.  Similarly,  if  /  only  appears  in  the  atomic  formulas  /  <  0  or  /  ^  0, 
the  point  x  =  a  does  not  need  to  be  tested. 

In  this  section  we  restate  Weispfenning’s  original  virtual  term  substitution  method, 
then  provide  a  different  substitution  for  infinitesimal  expressions,  one  that  uses 
fewer  atomic  formulas  and,  more  importantly,  fewer  distinct  polynomials.  In 
particular,  the  polynomials  in  the  resulting  formula  are  the  same  regardless  of 
whether  or  not  substitution  of  infinitesimals  is  required. 


3.1  Virtual  term  substitution  for  formulas  with  x-degree 
at  most  2 

Evaluating  F  at  a  point  a,  where  a  is  the  root  of  the  left-hand  side  of  some 
atomic  formula,  has  already  been  addressed.  Evaluating  F  at  —  oo  is  straight¬ 
forward.  Therefore,  evaluating  F  at  a  +  e,  where  a  is  the  root  of  the  left-hand 
side  of  some  atomic  formula,  is  what  remains.  Weispfenning  accomplishes  this 
by  considering  the  derivatives  of  polynomials  appearing  in  F,  and  thus  reduces 
the  determination  of  the  sign  of  a  polynomial  g  at  a  +  e  to  the  determination  of 
the  sign  of  g  and  its  derivatives  at  a.  In  short,  the  infinitesimals  are  removed, 
but  at  the  cost  of  introducing  new  polynomials. 


Theorem  6  (Weispfenning)  Let  f  =  ax2  +bx  +  c  and  let  g  =  agx2  +  bgx  +  cg. 
Assuming  a  ^  0  A  b2  —  4 ac  <  0, 

1.  g{ap  +  e)  =  0  ag  =  0  A  bg  =  0  A  cg  =  0 

2.  g(ap  +  e)  0  •<=>  ag  =  0  V  bg  =  0  V  cg  =  0 

3.  g(ap+e)  <  0  g(ap)  <  0 \/g(ap)  =  0A (gx(ap)  <  0 Vgx(ap)  =  0Aag  <  0) 
4 ■  g(ap  +  e)  <  0  <=$■  g(ap  +  e)  =  0  V  g{ap  +  e)  <  0. 

Assuming  a  =  0  A  b  ^  0,  all  of  the  above  holds  simply  replacing  ap  with  —c/b. 
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Let  Fap+e  denote  the  formula  obtained  by  carrying  out  the  a  V  0  substitutions 
for  each  element  of  F.  Let  F_c/b+e  denote  the  formula  obtained  by  carrying 
out  the  a  =  0  substitutions  for  each  element  of  F.  Let  -F_oo  denote  the  formula 
obtained  by  replacing  5  =  0  with  /\"=0  gi  =  0,  g  ^  0  with  V"=o  5*  A  and 

g  <  0  with  \/r=  o  <  0  A  A"=*+i  =  0) ,  where  5  =  gnxn  H - h  50- 

Let  /i  =  dja:2  +  biX  +  Ci  be  the  polynomial  occurring  in  the  ?’th  atomic  formula, 
fi  <Ji  0.  Let  A  =  b?  -  AaiCi  and  let  <^±1  =  ±  yjb'j  -  4alci)/(2ai)-  Let  I 

and  J  be  the  sets  of  indices  i  such  that  cr;  is  =,  <  and  <,  A,  respectively. 


Theorem  7  (Weispfenning)  3x  [T]  is  equivalent  to 

Vie/  (°i  =  ^  A  bi  V  0  A  F_c./b.  VaA0ADj>0A  {Fai,+1  V 
V 


VieJ  (a*  =  0  A  bi  V  0  A  F_Ci/bi+e  V  a*  V  0  A  A  >  0  A  (Fai  +1+e  V  Fai  _1+e)) 
V 


Notice  that  if  there  are  no  strict  inequalities  there  is  no  need  to  evaluate  at  a 
point  defined  by  infinitesimals.  Evaluation  at  infinitesimals  can  be  undesirable 
because  substitutions  for  g(ap  +  e)  <  0  and  g{ap  +  e)  <  0  require  substitutions 
for  gx(otp  +  e)  <0,  which  means  more  atomic  formulas  in  the  substituted  ex¬ 
pressions,  and  more  distinct  polynomials.  That  the  expression  substituted  for 
gx{ap  +  e)  <  0  really  contains  additional  polynomials  can  be  checked  by  simple 
calculation. 

•  If  /  =  ax 2  +  bx  +  c  and  g  =  ux2  +  vx  +  w,  then  r  =  av  —  ub,  a*  = 
2a2w  —  2 auc  —  bav  +  ub2,  and  R  =  4a(it2c2  —  2 ucaw  +  a2w2  —  vubc  — 
vbaw  +  av2c  +  wub2). 

•  If  /  =  ux2  +  vx  +  w  and  g  =  ax 2  +  bx  +  c,  then  r  =  —  ( av  —  ub),  a*  = 
2u2c—2uaw  —  vub+av2,  and  R  =  Au(u2c2  —  2ucaw  +  a2w2  —  vubc  —  vbaw  + 
av2c  +  wub2). 

•  If  /  =  ax2  +  bx  +  c  and  g  =  2 ux  +  v,  then  r  =  2 u,  a*  =  2(av  —  ub),  and 
R  =  4a(4u2c  —  2  vub  +  av2) 


This  shows  clearly  that  the  substitution  for  g(ap  +  e)  <  0  contains  the  polyno¬ 
mial  4u2c  —  2 vub  +  av2  which  is  a  part  of  neither  the  substitutions  for  g(ap)a 0 
nor  f(/3p)a0,  where  (3P  is  a  root  of  g.  Thus,  generically,  the  quantifier-free  for¬ 
mula  produced  by  Theorem  7  contains  more  atomic  formulas  and  more  distinct 
polynomials  when  infinitesimals  are  required  than  when  they  are  not. 
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3.2  A  simpler  substitution  for  infinitesimals 


In  this  section  we  give  a  simpler  substitution  for  infinitesimals  —  one  that  uses 
fewer  atomic  formulas  but  which,  more  importantly,  uses  the  same  polynomials 
as  are  used  without  infinitesimals.  The  key  observation  (see  Theorem  3)  is  that 
a*  =  —rfx(—s/r),  so  that  if  R  =  0  ^  r,  the  signs  of  a*  and  r  give  the  sign  of  fx 
at  the  common  root  of  /  and  g. 

Let  f,r,s  and  a*  be  as  before,  but  let  g  =  agx 2  +  bgx  +  cg.  Let  prem(/,  g)  = 
rgx  +  sg,  and  note  that  rg  =  —r  and  sg  =  —  s,  since  prern (f,g)  =  — prem (g,  /). 
Let  b*  =rg  =  -r  and  let  a*  =  2 agsg  -  bgrg  =  -rggx{-sg/rg)  =  rgx(-s/r). 


Theorem  8  g{ap  +  e)  <  0  is  equivalent  to 

R  >  0  A  a* a5  <  0  V  pras  <  0  A  ( a*as  <  0  V  R  <  0) 

V 

R  =  0/\(r  =  0/\  pag  <  0  V  a*  pi'  <  0  A  ra*  <  0  V  a*pr  <  0  A  a*  =  0  A  ag  <  0) 

under  the  assumption  a  /  0  A  52  -  4ac  >  0  for  p  =  +1  and  a  ^  0  A  b2  —  4 ac  >  0 
for  p  =  —  1 . 


PROOF.  g(ap  +  e)  <  0  is  equivalent  to  g(ap)  <  0  V  g(ap)  =  0  A  (, gx(oip )  < 
0  V  gx{ap)  =  0  A  ag  <  0).  The  first  line  of  the  substitution  formula  from  the 
theorem  statement  is  equivalent  to  g(ap)  <  0  by  Theorem  4  and  Corollary  1,  so 
we  focus  on  the  g(ap)  =  0  case. 

g(ap)  =  0  is  equivalent  to  a*pr  <  0  A  R  =  0.  If  R  =  r  =  0,  /  and  g  have  the 
same  roots.  In  this  case,  g{ap  +  e)  <  0  if  and  only  if  the  roots  are  distinct  and 
the  sign  ag  is  opposite  of  p,  or  ap  is  a  double  root,  in  which  case  the  sign  of  ag 
must  be  negative:  i.e.  b2  —  4ac  >  0  A pag  <  0  V  b2  —  4 ac  =  0Aas  <0.  If  p  =  +1 
or  b 2  —  4ac  >  0  this  can  be  simplified  to  pag  <  0. 

If  R  =  0  ^  r  A  a*pr  <  0  then  /  and  g  have  the  single  common  root  ap  =  —s/r. 
In  this  case,  a*  =  rgx(—s/r)  =  rgx(ap).  Therefore,  rgx(ap )  <  0  =>■  g(ap+e)  <  0 
and  rgx(ap)  >  0  =>  g(ap  +  e)  >  0  .  If  rgx(ap)  =  0,  ap  is  a  double  root  of  g  and 
a  simple  root  of  /,  so  g(ap  +  e)  <  0  <;==>  ag  <  0. 

So,  g{ap)  =  0  A  g(ap  +  e)  <  0  is  equivalent  to 

f?  =  0Ar  =  0A pag  <  0  V  R 
R  =  0  A  a*pr  <0Ar^0AaJ 

assuming  a  /  0  A  i)2  -  4ac  >  0  when  p 

p  =  -1. 


=  0  A  a*pr  <  0  A  ra*  <  0  V  . 

=  0  A  ag  <  0  (1) 

=  +1  and  a  /  0  A  i2  -  4ac  >  0  when 


If  R  =  0  yf  r,  a*  =  a*  =  0  implies  fx{—s/r )  =  gx(—s/r )  =  0.  But  this  means 
/  and  g  share  a  double  root,  which  contradicts  r  ^  0.  Thus,  R  =  0  A  a*pr  < 
0  A  r  /  0  A  o*  =  0  A  ag  <  Ois  false  when  a*  =  0,  and  we  may  simplify  it  to 
R  =  0  A  a*pr  <  0  A  a*  =  0  A  ag  <  0.  So  (1)  simplifies  to 

i?  =  0  A  (r  =  0A  pag  <  0  V  a*pr  <  0  A  ra*  <  0  V  a*pr  <  0  A  a*  =  0  A  ag  <  0)  . 

□ 


The  previous  theorem  showed  that  when  both  /  and  g  have  degree  two  g(ap  + 
e)  a  0  and  f{Pq  +  e)  p  0  can  be  rewritten  without  radicals  or  infinitesimals  using 
the  same  polynomials  that  are  used  in  rewriting  g{ap)  a  0  and  f(/3q)  p  0  without 
radicals.  The  next  theorem  shows  the  same  thing  when  /  has  degree  two  and  g 
has  degree  one. 

Theorem  9  If  g  =  ux  +  v,  where  u  yf  0,  and  let  j3  =  —v/u,  so  that  f(/3)  =  0. 
Let  R  =  R/(4a),  noting  that  the  division  is  exact. 

f(/3  +  e)  <  0  •<=>■  i?<0Vi?  =  0A  (— ua *  <  0  V  a*  =  0  A  a  <  0). 

Note  that  we  do  not  assume  that  a  ^  0. 

PROOF.  First  note  that  R/(4a)  —  res x{f,g)  =  res x(g,f)  =  u2f(f3),  so  that  the 
sign  of  R  is  the  sign  of  f(/3).  Then  note  that  prem(g,  /)  =  g  =  ux  +  v,  so  r  =  u 
and  s  =  v.  Thus,  by  Theorem  3,  a*  =  —rfx(—s/r)  =  — ufx(/3 ).  Since  a  /  0, 
sgn (fx(/3))  =  sgn (—ua*).  So,  /(/3  +  e)  <  0  is  equivalent  to 

J?<0VJZ  =  0A  ~ua*  <  0  V  R  =  0  A  a*  =  0  A  a  <  0  . 

/(/ 3)<0  /(/3)=0A/,(«<0  /  (P)=fx(P)—0Aa<0 


□ 


The  important  thing  about  this  theorem  is  that  it  shows  that  rewriting  /(/3+e)  < 
0  (or  <  0)  does  not  require  any  polynomials  that  are  not  already  required  in 
rewriting  g{ap)  a  0.  Notice  that  g(ap  +  e)  can  be  rewritten  by  specializing 
Theorem  8  setting  ag  =  0,  bg  =  u,  cg  =  v. 

Finally,  note  that  substitution  in  the  case  where  both  /  and  g  are  linear  is 
straightforward.  If  /  =  ax  +  b,  a  =  —b/a,  and  g  —  ux  +  v,  then  g(a)  a  0  is 
equivalent  to  a  res x(f,g)  &  0,  and  g(a  +  e)  <  0  is  equivalent  to  g{a)  <  0 V g(a)  = 
0A  u  <  0.  Thus,  for  any  combination  of  degrees  of  /  and  g  we  have  substitutions 
for  roots,  possibly  with  infinitesimals,  that  involve  at  most  the  coefficients  of  / 
and  g ,  discriminants,  pairwise  resultants,  first  principal  subresultant  coefficients, 
and  a*  and  a*.  Using  the  substitutions  from  Theorem’s  8  and  9,  we  must  use 
the  following  slight  modification  of  Theorem  7: 
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Theorem  10  3x[F]  is  equivalent  to 


Vie/  (ai  =  0  A  bi  ^  0  A  F_Ci/bi  V  a,;  ^  0  A  A  >  0  A  {Dai,+1  V  Fai  _±)) 
V 


V,;e./  | 

(  1 

at  =  0  A  bi  ^  0  A  F_Ci/bi+e  V  ai  ^  0  A 

f  Di  >0AFa.+1+e  \ 
V 

|) 

V 

V  ' 

V  Di  >  0A  Faii_1+e  ) 

7 

The  difference  between  Theorem  10  and  Theorem  7  is  that  instead  of  assuming 
Di  >  0  for  both  Fai  +1+e  and  Fai  _1+e,  we  assume  Di  >  0  for  Fai  +1+t^  and 
assume  Di  >  0  for  Fai  _1+e.  We  do  this,  of  course,  to  meet  the  requirements 
of  Theorem  8.  However,  it  makes  a  certain  sense,  because  now  the  Di  =  0  case 
is  covered  by  just  one  subformula.  Another  tangible  benefit  of  this  comes  from 
the  nice  way  we  can  substitute  for  f(ap  +  e)  <  0,  i.e.  /  evaluated  to  the  right 
of  one  of  its  own  roots.  This  can  now  be  rewritten  as  a  <  0  when  p  =  +1  and 
a  >  0  for  p  =  —  1,  because  the  possibility  of  a  double- root  when  p  =  —  1  has 
been  eliminated. 


4  A  different  view  of  virtual  term  substitution 


The  fundamental  question  to  be  addressed  in  the  quadratic  case  of  virtual  sub¬ 
stitution  is  this:  What  is  the  sign  of  g  at  root  ap  of  /  =  ax2  +  bx+c?  (Assuming, 
of  course,  that  a  /  0  and  b 2  —  4ac  >  0.)  The  answer  to  this  question  has  to 
be  expressed  as  a  combination  of  polynomial  equalities  and  inequalities  in  the 
remaining  variables  —  i.e.  without  x.  In  this  section  we  give  a  geometric  view 
of  how  this  is  done  and,  based  on  that  view,  suggest  alternatives  to  the  substi¬ 
tutions  given  in  Theorem  4. 


4.1  A  geometric  view  of  evaluation  at  roots  of  / 

Recall  that  R  =  4anresx(f,  g),  so  that  R  has  the  same  sign  as  the  product 
of  g  evaluated  at  the  two  roots  of  /.  There  is  a  geometry  to  this  problem  of 
determining  the  sign  of  g  at  ap  that  can  be  seen  quite  clearly  by  considering  r 
and  s  in  R  =  4 a(as2  —  srb  +  r2c)  as  variables  and  a,  6,  and  c  as  constant.  R  is 
the  product  of  two  lines  through  the  origin: 

~Cr  +  s)  (^f3Hr  +  s) 

For  a  specific  g:  If  (r,  s )  falls  on  the  first  line  g(a+)  =  0.  If  (r,  s)  falls  below  the 
first  line  g{a+)  <  0  .  If  (r,  s)  falls  on  the  second  line  g(a_)  =  0.  If  (r,  s)  falls 


10 


below  the  second  line,  g(a~)  <  0.  In  other  words,  the  sign  of  g(ap)  is  determined 
by  where  (r,  s)  falls  with  respect  to  these  two  lines.  Figure  4.1  illustrates  this 
for  a  specific  /. 


Figure  1:  The  region  in  which  g(a+)  <  0,  for  /  =  x2  +  x  —  1.  Points  in  the 
(r,  s)-plane  corresponding  to  two  different  g’s  are  shown. 


Of  course  we  do  not  want  to  evaluate  ~b+v)bg2 ~4—  r  +  s  and  ~b~^~4— r  +  s 
directly,  because  they  involve  radicals.  Instead  we  evaluate  R ,  a  multiple  of 
their  product.  However,  there  are  9  possible  combinations  of  sign  for  the  two 
linear  factors,  and  only  3  possible  signs  for  R.  Thus,  other  polynomials  need  to 
be  introduced  to  distinguish  between  regions  in  which  the  sign  of  R  is  the  same, 
but  the  signs  of  the  linear  factors  are  different.  This  is  precisely  the  role  of  r 
and  a*. 

Since  r  =  0  =>  R  >  0,  r  always  separates  the  two  regions  in  which  R  <  0. 
Since  4a2,  the  leading  coefficient  of  R  as  a  polynomial  in  s,  is  always  positive, 
dR/ds  =  4a(2as  —  rb)  =  4aa*  always  separates  the  two  regions  in  which  R  >  0. 
Geometrically,  r  and  a*  do  nothing  more  than  distinguish  between  disconnected 
regions  in  which  R  has  the  same  sign.  (Except  at  the  origin,  where  they  partition 
R  =  0  into  five  distinct  regions.) 

One  might  consider  whether  different  polynomials  could  be  used  to  distinguish 
these  regions.  This  would  provide  different  substitutions  than  those  from  The¬ 
orem  2.1  of  [6].  Since  the  original  substitutions  are  based  on  r  and  dR/ds ,  in 
the  following  section  we  examine  cases  in  which  s  and  dR/dr  can  be  used  as 
separating  polynomials  instead. 


4.2  Alternate  substitutions 

In  this  section  we  at  a  few  alternatives  to  the  substitutions  based  on  r  and 
a*  given  by  Weispfenning.  Whether  or  not  these  substitutions  may  be  applied 
depends  on  the  coefficients  of  /,  but  is  independent  of  g. 
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If  ac  <  0  we  note  that  R  is  the  product  of  two  lines  with  opposite  slopes  (see 
the  left  plot  from  Figure  4.2).  In  this  case  s  =  0  separates  the  R  <  0  regions 
just  as  does  2as  —  rb  =  0,  and  dR/dr  =  4a(2cr  —  6s)  =  0  separates  the  R  <  0 
regions  just  as  does  r  =  0. 


Figure  2:  Plots  of  R,  s  and  dR/dr  in  (r,  s)-space  for  ac  <  0  and  ac  >  0. 

If  ac  >  0,  R  is  the  product  of  two  lines,  both  of  which  have  positive  slope  if  ab  > 
0,  and  negative  slope  if  ab  <  0  (see  the  right  plot  from  Figure  4.2).  In  this  case 
s  =  0  separates  the  R  >  0  regions  just  as  does  r,  and  dR/dr  =  4a(2cr  —  6s)  =  0 
separates  the  R  <  0  regions  just  as  does  2 as  —  rb  =  0. 

Based  on  these  observations,  a  variety  of  alternate  substitutions  can  be  formu¬ 
lated  whose  applicability  is  dependent  on  the  signs  of  ac  and  R.  Figure  4.2  lists 
some  of  them  (note  that  c*  is  used  to  refer  to  2 cr  —  6s,  so  that  dR/dr  =  4oc*). 
Each  entry  has  been  verified  using  quantifier  elimination  Mathematica,  Red- 
log  and  Qepcad  B  all  verify  them  almost  instantly.  At  first  glance,  replacing  a* 
or  r  in  the  substitutions  from  Theorem  4  seems  to  require  assumptions  about 
the  sign  of  R  as  well  as  ac.  With  one  exception,  however,  a*  and  r  only  appear 
in  conjunction  with  the  required  sign  condition  on  R ,  so  that  really  only  the 
sign  of  ac  constrains  our  use  of  alternatives  for  a*  or  r.  The  one  exception  is 


If  ac  >  0  A  62  —  4ac  >  0  then 


If  ac  <  0  A  62  —  4ac  >  0  then 


If  ac  0  A  62  —  4ac  >  0  then 


R 

< 

R 

> 

R 

= 

R 

< 

R 

> 

R 

< 

R 

> 

R 

= 

R 

< 

R 

> 

= 

0 

0  =>  sgn(r)  =  sgn(a6s) 

0  =>  sgn(a*)  =  sgn(— a6c*) 
0  =>  sgn(a*r)  =  sgn(— c*s) 
0  =>  (r  =  0  s  =  0) 

0  =>  (a*  =  0  c*  =  0) 

0  ==>  sgn(r)  =  sgn(— ac*) 

0  ==>  sgn(a*)  =  sgn(as) 

0  ==>  sgn(a*r)  =  sgn(— c*s) 
0  =>  (r  =  0  <=>  c*  =  0) 

0  =>  (a*  =  0  <*=>  s  =  0) 
==>  sgn(a*r)  =  sgn(— c*s) 


Figure  3:  Alternate  substitutions. 
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substitution  (3)  of  Theorem  4,  in  which  pra6  <  OA a*  a6  is  not  explicitly  guarded 
by  any  sign  condition  on  R.  The  following  theorem,  whose  simple  proof  we  omit, 
states  that  R  =  0  is  actually  implicit  in  this  case,  and  therefore  that  we  may 
freely  replace  a*  and  r  in  the  substitutions  from  Theorem  4  with  the  alternatives 
given  in  Figure  4.2  based  solely  on  the  sign  of  ac. 

Theorem  11  If  R  >  0  implies  X  a*  a6  <  0  and  R  <  0  implies  Y 
pra6  <  0,  then  X  can  he  used  interchangeably  with  a*  a6  <  0  and  Y  can  be  used 
interchangeably  with  pra6  <  0  in  substitution  (3)  of  Theorem  f. 

In  asking  whether  these  substitutions  are  useful,  it  is  helps  to  consider  what 
happens  generically.  For  example,  when  /  =  ax 2  +  bx  +  c  and  g  =  ux 2  +vx  +  w, 
we  have  r  =  av  —  ub,  s  =  aw  —  uc,  a*  =  2 a2w  —  2 auc  —  bav  +  ub 2  and  c*  = 
2cav  —  cub  —  baw.  Clearly  in  this  generic  case,  both  s  and  c*  are  ’’better” 
substitutions  than  a*.  In  the  non-generic  case,  when  coefficients  are  constants 
or  are  algebraicly  related,  any  one  of  these  can  be  good  or  bad  substitutions. 
What’s  interesting  is  that  can  generate  each  of  them  and  choose  the  substitution 
that  works  best  for  each  /,  g  combination  in  the  context  of  the  problem  to  be 
solved.  Section  5.1  provides  an  interesting  application  of  this  approach. 


5  Examples 


This  section  steps  through  two  example  computations  —  the  first  involving  a 
quadratic  constraint,  the  second  involving  infinitesimals.  The  results  of  using 
the  original  substitutions  are  compared  with  using  the  improved  substitution  for 
infinitesimals  from  Section  3.2  and  the  alternate  substitutions  from  Section  4.2. 
One  difficulty  in  going  through  examples  of  virtual  term  substitution  in  detail 
is  that  the  formulas  are  so  large  that  they  are  hard  to  look  at.  We  will  endeavor 
to  ameliorate  this  by  showing  only  key  parts  of  the  substituted  formulas,  and  by 
performing  some  reasonable  simplifications  before  substituting.  Also,  we  note 
that  R  =  anresx(f,g),  so  where  R  appears  in  formulas  we  will  use  asR ,  where 
R  =  res  x{f,g). 


5.1  An  example  from  epidemiology 

Andreas  Weber  and  his  colleagues  have  been  working  on  applying  symbolic  tools 
to  investigations  of  epidemiological  models,  this  example  comes  from  his  work. 
In  considering  the  existence  of  an  ’’endemic  equilibrium”  for  the  SEIT  model 
[3],  a  system  of  ODEs  used  to  model  tuberculosis  and  other  diseases,  one  arrives 
after  straightforward  calculations  at  the  following  formula: 

3S  [.f(S)  =  0  A  —S  <  0  A  5  -  1  <  0] 
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where  /  =  uPi(p2  —  Pi)S2  +  (dPir2  —  d2p2  +  d2  Pi  +  Pirir2  —  dvfi 2  +  vPiqr2  — 
df32T2-\-dvf3\  —  P\vp2+  Pir\d)S + p2d{d-\-v +r2) ,  all  parameters  are  positive,  and 
/3i  >  /32.  Note  that  the  assumptions  on  the  parameters  imply  that  the  coefficient 
of  S 2  is  negative  and  the  coefficient  of  S'0  is  positive.  We  will  apply  virtual  term 
substitution  to  this  problem.  For  the  sake  of  brevity,  however,  we  will  not  give 
the  a  =  0  substitution  or  the  a+i  substitution,  both  of  which  produce  obviously 
unsatisfiable  subformulas.  Thus,  the  quantified  input  formula  is  equivalent  to: 

a^0Ad>0A  <7i(a_i)  <  0  A  g2(a_i)  <  0, 

where  g±  =  —S  and  g2  =  S  —  1.  Since  the  quadratic  and  constant  coefficients 
have  opposite  signs,  ay^0A.D>0is  always  true,  so  we  will  proceed  with 
gi{a-\)  <  0  A  g2(a_i)  <  0.  Since  p  =  —  1,  g\  and  g2  have  degree  1,  and  a 
is  always  negative  we  will  reduce  a*a s  <  0  A  asR  >  0  V  pra5  <  0  A  ( a* a 5  < 
0  V  asR  <  0)  to 

-a*<0A-i?>0Vr<0A  (-a*  <  0  V  -R  <  0) 

Following  Weispfenning’s  original  substitutions  restated  in  Theorem  4,  we  get: 

—  (d/?ir2  -  p2d2  +  d2(3i  +  P\r\r2  —  v(32d  +  iypiqr2  —  p2dr2  +  dvf3\  —  P\vP2 

+ Pirid)  <  0  A  -(dp2(d  +  v  +  r2))  >  0  V  —  1  <  0  A  [— (d/?ir2  —  P2d2 
+d2(3i  +  PiriV2  -  vfad  +  vpi qr2  -  P '2dr2  +  dvpi  -  Pivp2  +  Pirid)  <  0 
V  —  (dp2{d  +  v  +  r2))  <  0] 

A 

—  (2vf32  —  —  dPir2  +  P  2d2  —  d2P\  —  P\r\r2  +  vp2d  —  vPiqr2  +  P2dr2 

—dvpi  —  Pirid)  <  0  A  —(Pi(—PiV  +  r2d  +  d2  +  nr2  +  vqr2  +  dv 
+dr\))  >  0  V  1  <  0  A  (—(2vP2  —  PivP2  —  d/3ir2  +  p2d2  —  d2Pi  — 

Pirir2  +  vp2d  —  vPiqr2  +  P2dr2  —  dvpi  -  Pirid)  <  0  V  —  {Pi{~Piv 
+r2d  +  d2  +  rir2  +  vqr2  +  dv  +  dri))  <  0) 

Noting  that  dp2(d  +  v  +  r2)  is  always  positive  and  simplifying  away  inequalities 
involving  only  constants,  we  get: 

2  vPl  —  P1VP2  —  dP\r2  +  P2d2  —  d2P\  —  Pirir2  +  vp2d  —  vP\qr2  +  P2dr2 
—dvpi  —  P\r\d  >  0  A  —P\v  +  r2d  +  d2  +  r\r2  +  vqr2  +  dv  +  dr\  <  0 

However,  we  are  in  the  ac  <  0  case,  so  we  may  replace  a*  with  as  according 
to  Figure  4.2.  Since  si  =  0  and  s2  =  —  1  and  a  is  known  to  be  negative,  this 
alternate  substitution  is  well  worth  taking.  With  it  we  get: 

( — 0)  <  0  A  —  (dp2(d  +  v  +  r2))  >  0V 

—  1  <  0  A  [—(0)  <  0  V  —  (dp2{d  +  v  +  r2))  <  0] 

A 

-(-(-1))  <  0  A  —(Pi(—Pii/  +  r2d  +  d2  +  rir2  +  vqr2  +  dv  +  drq))  >  0 

VI  <  0A 

[—(—(—!))  <  0  V  —(Pi(—Piv  +  r2d  +  d2  +  rir2  +  vqr2  +  dv  +  drq))  <  0] 
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After  making  the  obvious  simplifications  we  get: 

—fd\v  +  r2d  +  d2  +  r\r2  +  vqr2  +  dv  +  dr\  <0  (3) 

The  final  simplification  of  (2)  to  (3)  is  not  trivial.  The  assumptions  on  the 
parameters  do  not  imply  the  positivity  of  the  extraneous  polynomial.  It  is  only 
those  assumptions  in  conjunction  with  —  p\u+r2d+d2  +rir2+vqr2+di,+dri  <  0 
that  imply  it.  This  is  a  simplification  that  Redlog’s  simplifier,  for  example,  is 
not  able  to  make. 


5.2  Substituting  infinitesimals 

Let  /  =  ax2+bx+ 1  and  g  =  ux2+vx—l.  We  consider  the  formula  3 x[f  <  0A g  < 
0]  under  the  assumption  a,  u  >  0.  First  we  will  follow  the  original  method,  then 
we  will  apply  the  improved  substitutions  for  infinitesimals  as  well  as  alternate 
substitutions  from  the  previous  section.  Rather  than  write  out  the  entire  formula 
here,  we  simply  write  out  the  set  of  polynomials  appearing  in  the  formula,  and 
show  one  representative  subformula,  the  substitution  for  g{ct-i  +  e)  <0.  In 
addition  to  the  coefficients  of  f  and  g ,  the  following  polynomials  appear  in  the 
formula  produced  by  Theorem  7 : 

Df  =  b2  —  4a,  Dg  =  v2  +  4 u,  r  =  av  —  ub ,  s  =  —a  —  u 
R  =  u2  +  2  au  +  a2  —  vub  +  bav  +  av2  —  ub2 

a*  =  —2a2  —  2 au  —  bav  +  ub 2,  a*  =  2 u2  +  2 au  —  vub  +  av2  (4) 

R9x  =  4u2  —  2 vub  +  av2,  Rfx  =  —4a2  —  2bav  +  ub2 
c*  =  2 av  —  ub  +  ab,  c*  =  av  —  2 ub  —  vu 

The  original  substitution  for  g(a-i  +  e)  <  0  is: 

gx(a-i)<0 

/  2 ra  <  0  A  aR9x  >  0 

I  V  —  2 ua  <  0A 
\  (2 ra  <  0  V  aR9x  <  0)) 

V 

— 4itr  <  0  A  aRgx  =  0  A  2 u 

V 

gx(a_i)=0A2u<0 

Although  the  other  substitutions  are  not  shown,  it  should  be  clear  that  the 
entire  formula  is  constructed  out  of  the  polynomials  in  (4)  and  a,  b ,  u,  v. 

The  substitution  given  in  Section  3.2  gives  the  following  for  g(a- 1  +  e)  <  0: 

(a*  <  0  A  H  >  0  V  -r  <  0  A  (a*  <  0  V  ~R  <  0)) 

V 

i?  =  0A(r  =  0A  —u  <  0  V  -a*r  <  0  A  ra*  <  0  V  -a*r  <OAa*=0A«<0) 


g(a_i)<0 


a*  <  0  A  R  >  0 
V  —  r  <  0A 
(a*  <0VR<0)) 


g(a~i)=0 
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Although  the  other  substitutions  are  not  shown,  it  should  be  clear  that  the 
entire  formula  is  constructed  out  of  a,  b ,  u,  v  and  the  polynomials  in  (4)  minus 
Rgx  and  Rfx. 

Clearly  /  falls  in  the  ac  >  0  case  and  g  falls  in  the  ac  <  0  case  discussed  in 
Section  4.2.  Thus,  we  may  choose  alternate  substitutions  given  there.  This 
leads  to  a  substitution  for  g(a- 1  +  e)  <  0  of: 

(— abc *  <  0  A  R  >  0  V  — abs  <  0  A  (— abc *  <  0  V  R  <  0)) 

V 

R  =  0  A  (s  =  0  A  —u  <  0  V  sc*  <  0  A  sc*  <  0  V  sc*  <  0  A  us  =  0  A  u  <  0) 

Whether  or  not  this  is  ”  better”  than  the  previous  formula  depends  on  what 
subsequent  computation  is  desired.  It  is  interesting,  however,  that  it  is  trivial 
to  deduce  that  the  input  assumptions  a,  u  >  0  implies  s  >  0,  which  then  con¬ 
siderably  simplifies  the  formula.  It  is  also  interesting  that  after  removing  all 
polynomials  that,  by  inspection,  never  vanish  given  a,  u  >  0,  this  final  version 
contains  only  2  polynomials  that  are  not  linear  —  R  and  Df.  In  contrast,  the 
original  contains  4  non-linear  polynomials.  The  potential  advantage  to  alter¬ 
nate  substitutions  is  that  a  program  may  quickly  examine  the  alternatives  and 
decide  whether,  as  in  this  case,  one  offers  advantages  over  the  other. 


6  Improved  bound  for  virtual  term  substitution 


The  fact  that  res x{f,g)  =  (a*2  —  b*2D)/(Aan)  allows  one  to  tighten  the  most 
general  degree  bound  given  in  Corollary  2.2  of  [6].  Suppose  that  M  is  the 
maximum  total  degree  of  any  polynomial  in  the  input,  and  that  d  is  the  greatest 
degree  in  x  of  any  polynomial  in  the  input. 

Theorem  12  The  highest  degree  of  any  irreducible  factor  of  a  polynomial  ap¬ 
pearing  in  the  formula  produced  by  Theorem  2.1  of  [6]  is  (d  +  2 )M  —  2d. 

PROOF.  Note  that  the  total  degree  of  the  coefficient  of  xm  is  at  most  M—m.  The 
candidates  for  the  highest  degree  factors  are  b2  —  4ac,  b*,  a*  and  a*2  —  b*2c.  The 
degree  of  b2  —  4 ac  is  clearly  bounded  by  2 M  —  2.  a*2  —  b*2c  is  the  determinant 
of  the  Sylvester  matrix  for  /  and  g ,  and  r  and  s  are  given  by  minors  of  the 
Sylvester  matrix.  We  will  show  explicitly  that  the  largest  irreducible  factor  of 
a*2  —b*2c  has  total  degree  at  most  (d+2)M  —  2d.  A  similar  approach  shows  that 
r  and  s  have  total  degrees  at  most  dM  —  2d  +  1  and  dM  —  2d +  2,  respectively. 
Thus,  a*  =  2as  —  br  lias  total  degree  at  most  max((M  —  1)  +  dM  —  2d+  1,  (M  — 
2)  +  dM  -  2d  +  2)  =  (d  +  1)M  -  2d. 

res x{f,g)  =  (a*2  —  b*2D)/(4an),  the  degree  of  the  largest  irreducible  factor  of 
a*2  —  b*2D  is  bounded  from  above  by  the  degree  of  res^f?,  /).  Assume  g  has  the 
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maximal  x-degree  d.  The  rows  of  the  Sylvester  matrix  for  g  and  /  correspond  to 
xg,  g,  xn~2f,  xn_1/, . . . ,  x°f.  The  determinant  is  the  sum  of  all  products  of  one 
element  from  each  row  and  each  column.  Consider  choosing  elements  to  form 
such  a  product.  Suppose  that  i  and  j,  i  <  j,  are  the  indices  of  the  entries  chosen 
from  the  first  two  rows.  From  columns  1 , ,i  —  1  we  must  choose  the  a  entry 
in  order  to  get  a  non-zero  product  (a  has  degree  at  most  M  —  2).  From  columns 
j  +  1, . . . ,  d  +  2  we  must  choose  the  c  entry  to  get  a  non-zero  product  (degree 
M).  The  submatrix  remaining  after  all  these  choices  is  tridiagonal  with  a's 
below,  b's  on  and  c’s  above  the  diagonal.  Any  entry  chosen  above  the  diagonal 
must  be  matched  with  an  entry  below  the  diagonal,  so  the  average  total  degree 
is  (M  —  1).  The  product  of  the  two  entries  from  the  first  two  rows  has  degree 
2 M  —  2d  +  i  +  j  —  3.  Thus,  any  term  in  the  determinant  has  degree  at  most 

(i—l)(M—2)+(j—i—l)(M—l)+(d+2—j)M+(2M—2d+i+j—3)  =  (d+2)M-2d. 

□ 


Corollary  2.2  of  [6]  gives  a  bound  of  (2d+2)M  —  2d  on  any  polynomial  appearing 
in  the  formula.  The  new  bound  is  approximately  a  factor  of  two  improvement, 
although  of  course  it  is  a  bound  on  the  size  of  irreducible  factors.  The  bound 
from  Corollary  2.2  also  assumes  that  the  total  degree  of  /  is  not  more  than  the 
maximum  total  degree  of  p{F).  The  above  analysis  makes  no  such  assumption. 


7  Conclusion 


This  paper  provides  an  analysis  of  the  polynomials  appearing  in  Weispfenning’s 
method  of  quantifier  elimination  by  virtual  term  substitution.  Based  on  this 
analysis,  and  simpler  substitution  is  given  for  the  evaluation  of  a  formula  at 
x  =  a  +  e,  where  a  is  the  root  of  a  quadratic  polynomial  and  e  is  a  positive 
infinitesimal.  The  paper  proceeds  with  a  new  view  on  why  certain  polynomials 
appear  in  substitutions  and,  based  on  this,  proposes  alternate  substitutions. 
These  alternatives  are  not  always  applicable  but,  when  they  are,  they  allow 
for  an  implementation  of  virtual  term  substitution  that  can  choose  amongst 
alternatives  in  order  to  produce  simpler  formulas.  Both  of  these  improvements 
are  aimed  at  helping  reduce  the  complexity  of  the  result  of  quantifier  elimination 
by  virtual  term  substitution,  which  is  the  method’s  biggest  problem. 
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